Derealization in an open one-dimensional chain in an imaginary vector potential 
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We present first results for the transmittance, T, through a ID disordered system with an imag- 
inary vector potential, ih, which provide a new analytical criterion for a derealization transition in 
the model. It turns out that the position of the critical curve on the complex energy plane (i.e. the 
curve where an exponential decay of (T) is changed by a power-law one) is different from that ob- 
tained previously from the complex energy spectra. Corresponding curves for (T n ) or (lnT) are also 
different. This happens because of different scales of the exponential decay of one-particle Green's 
functions (GF) defining the spectra and many-particle GF governing transport characteristics, and 
reflects higher-order correlations in localized eigenstates of the non-Hermitian model. 
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Non-Hermitian models with disorder have attracted a 
great deal of attention, e.g. in context of open quantum- 
mechanical (or optical) systems connected to reservoirs 
or random walks in disordered media Q , or Gaus- 
sian ensembles of non-Hermitian random matrices 
The discovery by Hatano and Nelson j|] of a dereal- 
ization transition in a simple Anderson model with an 
imaginary vector potential has led to a new burst of ac- 
tivity in this area j?J and also focused research interest 
on this particular model [p[-pl| and related problems |L2| . 

The numerical results on the complex spectra of the 
model [p| have been confirmed analytically, both for the 
zero- D [p| and ID || cases. The fact that the transi- 
tion into complex spectra in this model is indeed a de- 
localization transition has been convincingly illustrated 
by direct numerical investigations Jll| of the product of 
left and right eigenfunctions corresponding to the same 
complex eigenvalue. This product, which should enter 
into any physical observable, has shown to be delocalized 
JTT| |, in spite of the fact that eigenfunctions by them- 
selves still look localized even in complex spectra [ p0| . 
It demonstrates that the description of the localization- 
delocalization transition may depend on the choice of the 
localization criterion applied, in contrast to the standard 
Hermitian model. This makes it interesting to look at 
some analytical criterion of the transition. 

One of the best ways to address this problem is to cal- 
culate the transmittance T or, equivalently, the Landauer 
conductance. Short of a direct analytical calculation of 
the wavefunction distribution, an analytical dependence 
of T on the system size L provides the most straightfor- 
ward criterion of the derealization transition. It allows 
one to investigate how the localization length £ depends 
on the particle energy and how it diverges in approaching 
the transition. It is well known that the transmittance 
through a ID Hermitian disordered wire decays as 



T = a(L) exp 



(i) 



indicating that all the wavefunctions are localized with 
the localization length £ (depending on energy). The 



existence of the derealization transition should be man- 
ifested by the divergence in £ at the mobility edge. 

In this Letter, we present the analytical results for the 
transmittance T which confirm the existence of the der- 
ealization transition in the ID model introduced in [0. 
However, properties of the transition and even the po- 
sition of the mobility edge on the complex energy plane 
turn out to be different from those expected from the 
knowledge of the complex spectra alone. Before describ- 
ing this, note that T can only be analyzed for an open 
system. If one considered a system with free boundaries, 
the imaginary vector potential, h, could be wiped out by 
a gauge transformation. We consider a disordered sam- 
ple attached to ideal leads, with the logarithmic deriva- 
tive of the wave function being continuous at the bound- 
aries. Although h can be eliminated from the sample by 
the gauge transformation, the requirement that the wave 
functions remain finite in the leads imposes a nontrivial 
constraint on the class of allowed functions, similar to 
that imposed by the periodic boundary conditions. 

Although the open geometry for the model with an 
imaginary vector potential could not be directly related 
to the problem of depinning of flux lines in a supercon- 
ductor with columnar defects fl3f| , which was the original 
motivation for this model Q , the subsequent interest has 
arisen due to a general character of the observed hyper- 
sensitivity of a disordered system to the presence of even 
a small non-Hcrmiticity. Thus we believe that the ana- 
lytical description in the case of the open geometry can 
be very useful, even if not directly related to the original 
model of depinning [ p"3| . 

The model is described by the Schrodinger equation 



Hip{x) 



dx 



v(x) 



tl>(x) = zil>(x) . (2) 



Here ifj(x)=t/jfi(x) is the right eigenfunction of TL corre- 
sponding to the complex energy z\ in this model the left 
eigenfunction at the same energy is ipc{x; h) = ip^(x; —h) 
||. The random potential, v(x), is chosen to be Gaussian 



white-noise with zero mean and variance 



It vanishes 



1 



in the ideal leads attached to the sample at x = and 
x = —L. The imaginary vector potential, ih, is homoge- 
neous inside the sample and make take a different con- 
stant value (not necessarily 0) in the leads. The logarith- 
mic derivative of tp(x) is continuous at the boundaries. 

In the non-Hcrmitian case considered, we must find T 
as a function of complex energy z. Then the localization- 
delocalization transition would reveal itself via the diver- 
gence of £(z) at certain values of z (the mobility edge). 
Furthermore, for the values of z at which the wavefunc- 
tions are delocalized, the transmittance is expected to 
have a non-exponential dependence on L. 

However, solving the standard scattering problem for 
the Hamiltonian y2|) , we find that T always remains expo- 
nential as in Eq. (fil), albeit with £ increasing as a function 
of h. On the face of it, this result is quite surprising. It 
means that, were one able to have an experimental re- 
alization of a ID wire in the presence of the imaginary 
vector potential, measuring a current through it would 
not reveal an insulator-to-metal transition, at least as a 
function of the wire length. The reason for this is rather 
simple: the energies of incident plane waves on the com- 
plex z plane do not overlap with the energies of the de- 
localized states. Although these states exist, they make 
only an exponentially small contribution to the transmit- 
tance. Therefore, we need to generalize the scattering 
technique to be able to detect the delocalized states. 

We will define the scattering amplitudes, t and r, via 
the asymptotics of the wave function: 

^ ' ~ \ e lk +( x + L ) + re -ik-(x+L) X <-L. ^ ' 

We suppose, for a moment, that the imaginary vector po- 
tential h is the same both inside and outside the sample. 
Then the generalized wave- vector, k±, is defined by 

k± = \[z =p ih = k ± i(n — h) , 

with k and k being the real and imaginary parts of y/z 
(a branch with n > is chosen). Both the incident and 
transmitted waves are not divergent at x — > oo only pro- 
vided that k + is real, i.e. on the curve So in the z-planc, 

S : n = h 5Rcz= (3mz/2/i) 2 -/i 2 , (4) 

which happens to be also the DoS support curve in the 
absence of the impurity potential v. The DoS support 
inside the sample (v ^ 0) is, however, entirely different. 

To illustrate this, we first consider a simple case of the 
one-particle Green function (GF), G(z), which is the GF 
of the Schrodinger equation (|J). A straightforward cal- 
culation in the absence of the disorder, v(x) = 0, yields 

r> I n ^ 1 ( 6{n-h)e ik + x x>0 
^ 0[X ' U]Z) ~2i^ \e- ik - x -6{h- K )e ik + x x<0 W 

As a function of z, Go is discontinuous at k = ft which 
means that all the eigenvalues lie on the curve So in the 



z-plane defined by Eq. (|f|). This can easily be seen from 
the well known formula for the density of states (DoS), 
v{z) = 7r -1 <9^G(0, 0; z). Thus, the DoS support lies on 
the parabola along the 3?e z axis on the z plane. 

In the presence of the disorder it is straightforward 
to show that in the quasi-classical regime, k 3> K, £~ x , h, 
the ensemble-averaged GF of Eq. (||) remains of the same 
form as Go, Eq. (^), with the only change 

k± -> fc± ±i/2£(k), (6) 

the same as for the Hermitian problem, with £ = £(k) = 
2k 2 /it 2 being a mean free path. The discontinuity line of 
the GF defines the DoS support inside the sample: 

Si : k = max{/i- l/2£(k), 0} . (7) 

The DoS support lines inside and outside the sample 
(with the same imaginary vector potential, ih, in both 
regions), Eqs. (Q) and (Q), do not overlap at any disor- 
der. If we put h = outside the sample, then all states 
are on the line k = which only intersects Si at a sin- 
gle point. In order to detect delocalized states inside the 
sample, k + in the incident wave, Eq. (^) should be al- 
lowed to take on any complex value. To tune one 
could formally apply the imaginary potential iho outside 
of the sample. Then, by changing ho-, one would shift the 
DoS support line (Q) in such a way that it will intersect 
with Si (see Fig. 1). Note that any physical quantity, 
like conductance cx T, would be defined by the ratio of 
the wave functions on the sample boundaries, and thus 
will be totally unaffected by h . 

One would expect that Si also defines scattering prop- 
erties of the disordered system, in particular the trans- 
mittance T and thus the mobility edge where — > 
oo. Indeed, Eq. (|7]) defines a parabola on the z plane 
which is an exact counter-part of the bubble-like curve 
analytically found in for the lattice variant of the 
model. However, the situation turns out to be not that 
simple. What we have found is that the transmittance 
cumulants, (T n ), are characterized by the whole set of 
the critical curves, C n (see Fig. 1). 



il K/h 




h^(k)=l/16 =1/4 =1/2 k 

FIG. 1. Phase diagram in y^-plane: the DoS support lines 
for a pure system, So, and a disordered system, Si; the critical 
curves for (InT), S 2 , and (T), d; £(k) = 2k 2 /u 2 , Eq. (§). 
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This multitude of the critical curves is due to the ab- 
sence of self-averaging in ID It is InT which is 
normally-distributed in the Hcrmitian case. The lognor- 
mal distribution of T leads to all the moments T n to 
have the same decay exponent —L/£ (which is not scaled 
with n), where £ is the localization length. It is the two- 
particle GF which decays oc e~ L ^ and thus contains the 
information about the localization, while the decay of 
the one-particle GF, Eq. @, is only due to the loss of 
the momentum direction in elastic scattering, which is 
the same in any dimensionality. The self-averaging InT 
is governed by its own exponent different from that in 
the one-particle GF by substituting the transport scat- 
tering time, £ t i = 21, for £. In the presence of ih, the 
exponential decay is replaced by e^ n/t_1 /^ L for T n , or by 
e {h-i/2t)L for the on e-p a rticle GF and (h-l/2£ tI )L for 
InT. Then for h>l/£ (or 1/21, l/2£ tT , respectively) 
no small readjustment of localized states could prevent 
a mismatching of boundary conditions for a closed sys- 
tem (or an exponential increase in one of the leads for an 
open system). This defines a set of critical curves, Fig. 1, 
which arc different for all these quantities. 

Now we outline our procedure of rigorous analytical 
deriving of these results. By introducing vector VE", 



± 



1 dV> 
2iy/z dx 



(8) 



we reduce the Schrodinger equation (g) to the first-order 
matrix equation: 

— =/i$, fi(x) = h + i^fzaz + -—=(a 3 + ia 2 ) , (9) 

Q.X Z 

where Cj are the Pauli matrices. The corresponding 
transfer-matrix, m, is defined by the relation 

= m(x,x')^(x') (10) 

and obeys the following first order differential equation 

dm(x, x') 



dx 



n(x) ■ m(x, x') , 



(11) 



which is subject to the initial condition m(x,x) — 1. 

Taking into account that the values of ^ at the bound- 
aries, ^ T (-L) = (l,r) and * T (0) = (t, 0) are related 
via Eq. ( jlC)| ) by m(— L, 0), we express the reflection and 
transmission amplitudes in Eq. (j^) as follows 



t = 



mn (-£,0)' 



m 2 i (-T,0) 
m n {-L,0) 



(12) 



From Eq.(ll) we can extract two coupled equations for 
the scattering amplitudes r and t: 



dr „. r- v(—L) , ,2 
— = 2i^zr+^- 7 l 1 + r f, 
dL 2iy/z 

dt ,, v(—L) , . 

_ = ( fc + V5) t+ ^ t (l + r) 



(13) 



which obey the boundary conditions at the right end of 
the sample, r(0) = and t(0) = 1. By performing the 
ensemble averaging, one can easily derive equations for 
the mixed moments (T n R n ). As T + R ' ^ 1 in the non- 



Hcrmitian case, R : 



cannot be simply excluded. For 



our purposes, it is sufficient to consider only conditional 
averages, P n (L) = (T n (L, R)S (R(L) - R)) so that the 
moments (T n ) could be found by integrating over all R. 
As the reflectance R is an auxiliary quantity, it is conve- 
nient to introduce a new variable, x, by R = tanh 2 (x/2). 
In the absence of the imaginary potential, L/x is equal 
to the localization length [19] . The conditional probabil- 
ity P n (T,x) in new variables (with r = L/21) obeys the 
following Fokker-Planck equation, 



dPn 

dr 



d_ 

dx 



d 9fL 



dx dx 



(14) 



where the initial condition is P n (r = 0) = S(x), and 



2 X 



V n (x) = n tanh n — 4n(/t — h)£ , 

x ( 15 ) 

51„ (x) = 4:k£ cosh x — In sinh x + 4n In cosh — , 



which is derived from Eqs. (13) for k 3> k, h,£~ x in a 
way similar to that for the imaginary scalar potential . 
The moments of T can be found by the integration: 



(T») = / d X P n {x,r). 
Jo 



(16) 



As InT is well known to be a self-averaging quantity in 
ID, the Lyapunov exponent, 

A=-lim<^>, (lnT)=lim^l, (17) 

found via the standard 'replica' trick, gives the 'best' 
representation of the inverse localization length. 

Equation ([l4|) cannot be solved exactly and we map it 
onto the imaginary-time Schrodinger equation: 



dr 



dx 2 



+ U n (x)<S>n- 



(18) 



Here ® n {x, T ) = Pn{x, t) exp[ft n (x)], and the effective 
potential U n = (^) 2 /4 - Q'^/2 - V n reduces to 



Unix) 



4 4 sinh x 



-4nh£ 



+ A(k£) 2 sinh 2 x + in£{n ~ 1) coshx • 



(19) 



Although the eigenfunctions & n {x) of Eq. (Jig), where 
®n(Xi T ) = Ylj ^n(x) e~ £ " T , cannot be found exactly, the 
form of the effective potential, Eq. ( |l9| ) , makes possible 
to find the long-r limit of <& n (Xi T )i an( l thus of P n (Xi T )- 
A sharp increase in U n at x ~ ln(l/«^) makes the spec- 
trum of Eq. (18) discrete, with a gap of order 1 separating 



3 



the ground and excited states. Therefore, only the lowest 
eigenstates contribute to <E>„(x,t) and thus to (T n ) for 
t = L/21 3> 1. It is easy to verify that for n = the 
ground state has the energy e° = and the eigenfunction 



[/ d x e ~ 



•n (x) 



1/2- 



(20) 



For n -C 1, by treating the n-dependent part of the po- 
tential (|l9| ) as a perturbation, one finds $2 w \&q and 
e° = n[l + 4(k - h)£]. This leads, via Eq. (0), to exact 



1 £ U 

A = lim lim — = 2(k - h) 

2i T^oon^o n 



1 

27 



(21) 



The Lyapunov exponent X(z) vanishes (i.e. the localiza- 
tion length diverges) on the curve (see Fig. 1) 



l/U{k). 



(22) 



This curve, as expected, could be obtained from Si, 
Eq.(@), by substituting there £ t r = 21 for t. 

In contrast to (hiT), the moments (T n ) cannot be 
found exactly for n > 1. However, as rf < 1, approx- 
imate eigenfunctions can be found by separating % < 1 
and x ?i ln(/c£) scales and matching appropriate solu- 
tions. Thus we find the moments with n -C (k^) -1 : 

(T n (r; z)) = c„ (1/r + l/r c ) 3/2 e-[i- 2ll " + ^ r , (23) 

where r c (n) = 7r" 2 [ln(l/2/c£) - V>(ri- 1/2)] 2 , ip is the di- 
gamma function, c„ = [(2n-3)!!] 2 vr 5 / 2 2 1 - 2 ™ [(rc-l)!] -1 . 
For moments with n 3> (k^) _1 , the preexponential factor 
is proportional to \J n£/n, while the exponent becomes 
1/4 - 2n£(h - k). One finds from Eq. @ that the n-th 
exponent vanishes on the curve C n (which is drawn for 
n = 1 in Fig. 1 and is qualitatively the same for any n): 



2n£(k) 



3 -[2nM(fc)-l/4]" 1/2 



The corresponding transmittance moment on this curve 
close to the point where C n hits the &;-axis on the k-n 
plane, i.e. for h£(k) — l/8n <C 1, is given by 



T„(r; k) = c„ (r^ 1 + 2nM(k) - 1/4) 



3/2 



(25) 



Such a power- low behavior, (T n ) cx {IjLfl 2 , is typi- 
cal for a critical regime at the metal-insulator transition 
point. What is unusual is that each moment reaches the 
criticality at a different point in the z plane. 

The existence of a set of different critical curves means 
that the delocalization transition, say in (T), happens in 
the point in the z plane where the DoS support for a 
corresponding closed system, Si, is still on the real axis 
(Fig. 1) and therefore all products of the left and right 
eigenfunctions are still localized. This manifests the exis- 
tence of some higher-order correlations between the local- 
ized states, similar to the anti-correlations between the 



amplitudes of ipn(z) and ipc( z ) which ensures the delo- 
calization of their product on the curve Si, despite each 
of these function by itself is still 'localized', i.e. ooccupies 
only a very small part of the sample PJ. As the prop- 
erties of the localized states should not be strongly de- 
pendent on boundary conditions, these higher-order cor- 
relations between them could also exist for the closed 
disordered systems with the imaginary vector potential. 
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